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The connection of the Callen-Welton and Nyquist fluctuation-dissipation relations 
is considered for plasma-like classical and quantum systems. The conditions for ap- 
pearance of the dissipative parameters in the equilibrium current-current correlation 
function are investigated. The paper presents the arguments for the restrictions 
of the Nyquist theorem and against violation of the Callen-Welton theorem in the 
quantum case. 

I. INTRODUCTION 

The Callen-Welton theorem (CWT) is a general and rigorous result of statistical physics 
(see e.g. yfl). Some ideas have been suggested to show that the existing form of CWT 
in the quantum case is not correct. One of the arguments is based on a consideration of the 



oscillatory electric circuit. As is well known, the classical Nyquist formula 
electromotive force (emf) £ 2 has the form: 



for a random 



4 = 2TReZ(u), (1) 

where T is the temperature in energetic units and Z(u) is the impedance Z[iS) of a linear 
electric chain with the components arranged serially (for example). For a quasi-static field, 
changing slowly in time, Z(uj) is a function of the resistivity R, inductance L and capacity 
C for the electric RLC circuit (the speed of light is taken equal to unity c = 1): 
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In the general case the function Z[iS) is the Fourier- component of the linear operator Z, 
which connects the functions e(t) and J(t) in a linear relation: 

e = ZJ (3) 

or in Fourier-components: 

e w = Z(u)J u . (4) 
For the current fluctuations one derives: 

I J u T= 2T 2, (5) 

where R' and R" are the real and imaginary parts of the resistivity R{uj), respectively. Let 
us now calculate the full average square of the magnetic energy of the current LJ 2 /2: 

hi=l r 3 — i j \ 2 =lt r — - (6) 

2 2 7_ 00 2vr 1 w| J_ 00 27r RI 2 + ( R „ + ujL _J_y U 

For the case of ^-independent (and, therefore, real) resistivity R , which was considered 



in Q], Eq. (jHJ) can be written in the form: 



LJ 2 m f°° d£ 1 



in which ^ = Lu/Rq. The integral in Eq. (J7|) can be calculated analytically (or by use of, 
e.g., the mathematical program "Maple") and is equal to 1/2. There is no dependence of 
the integral (J2J) on the parameter LjCR\\ therefore it is also independent of the dissipative 
parameter R Q . The analytical calculation of the integral (JJJ) for the considered classical case 
of cu-independent active resistivity R is presented in the Appendix. 

II. FORMULATION OF THE PROBLEM, THE DISSIPATIVE AND PHYSICAL 

PARAMETERS 



DD 



The argumentation |2j, |3j against the quantum generalization of the Nyquist theorem 

el = 2RT u} - T^ = —coth— (8) 

is essentially based on the notion, that the magnetic energy of the current (see Eq. ©), 
as well as other correlation functions, calculated over the equilibrium distribution, cannot 
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include the dissipative parameter [2|,|3j]. More concretely, this statement has been formulated 
as the requirement for the integral of the magnetic energy LJ 2 /2 to be independent of the 
active resistivity Rq. This statement does not apply to the case of frequency-dependent 
Tw, when the magnetic energy is determined by an equation similar to Eq. (jHJ. In the 
quantum case, however, the integral in Q contains T w , instead of the temperature T in 

On 

the classical case, and depends on Rq. According to (2, |3J it is the reason to reconsider 
the quantum Nyquist formula and the other quantum fluctuation-dissipation relations or 
theorems (FDT). We want to stress that in P the independence of the charge fluctuations 
in the electric circuit has been considered on the basis of thermodynamic requirements only 
for the quasi-static case (hu) < T), e.g., in the classical case. 

At first we have to consider and analyze the requirement of independence of the equi- 
librium correlation functions of the "dissipative parameters", as the basic one for recon- 
sideration of the quantum FDT, in more detail. The simplest way to do that is to analyze 
microscopically the systems with Coulomb interaction, possessing a well-known Hamiltonian. 
It is evident, that the "dissipative parameters", such as, e.g., resistivity, themselves are func- 
tions of the fundamental non-dissipative parameters, namely, the mass m a , the charges e a , 
the particle density n a (for the finite samples the characteristic sizes can appear) and also 
the temperature T. Therefore it is necessary to clarify this statement by considering the 
possible dependence of the equilibrium correlation functions on the dimensionless parame- 
ters. We can suggest, for example, that for plasma-like systems the integral over frequencies 
of the current- current correlation function I = J 2 is of the form I = TQ 2 f(T, A, Z), where 
Z, T ~ e 2 n 1 / 3 /T and A ~ h 2 n 2 ^ 3 /mT are the plasma frequency, ion charge, interaction 
and degeneration parameters, respectively (the derivation of this statement is presented in 
Section III). 

Instead of T and A it is sometimes more convenient to use, as independent parameters, 
e.g., the parameters T and // = d/a® = T/A, where d = n~ x l 3 and ao is the Bohr radius. In 
a finite system additional dimensionless parameters connected with the sizes of the sample 
can be present. In some particular case this integral can be independent of T and has a form 



It seems natural to suggest that the statement on independence of the equilibrium corre- 
lation functions of the "dissipative parameters" in Q] can be formulated (in the case of the 



/ ~ m 2 p . 
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Otherwise, if there is dependence on the other parameters mentioned above (at first on 
T, for almost classical plasma) the dissipation manifests itself in the correlation function. 
This conclusion follows, in particular, from the formula for the collision frequency of the 
Lorentz plasma vq = 7rfi p r 3 / 2 L/8\/2 of a weakly non-ideal plasma with Z — 1. Here L is 
the Coulomb logarithm, which for a weakly non-ideal (r < 1) classical plasmas can be taken 
(with logarithmic accuracy) as L = ln(l/r). The frequency v is, of course, the " dissipative" 
parameter in the sense of ^. Therefore, it is natural to suggest that the dependence on 
"dissipative" parameters for J 2 means in fact the dependence on V. In any case, it is 
necessary to suggest which kind of dependence of the correlation functions, e.g. J 2 , on 
the selected parameters can, in principle, exist. Only on this basis can be proved that the 
FDT are rigorously correct or that in some cases, e.g., for quantum systems an incorrectness 
exists. 

We turn now to the Callen and Welton theorem (CWT) The CWT, strictly speaking, 
is formulated on the basis of statistical consideration of a system with a well- determined 
Hamiltonian and, therefore, can be applied to the Coulomb system with an arbitrary inter- 
action parameter. 

How to apply these results to the "macroscopic" electrodynamics, which are described 
in terms of macroscopical parameters as inductance and capacity? The answer is not quite 
clear, since an exact determination of correlation functions for these cases is lacking - as far 
as we know. 

It should be noted that the work (per unit time) produced by emf e under the moving 
charges in the linear circuit is equal to the sum of Joule heating, change of the magnetic 
field and the energy of electric field in the capacitor 

eJ = RJ i + -[— + — . (9) 



dt V 2 2C J 

This relation leads to the expression (J2J) for the frequency dependent impedance. At the 
same time the applicability of this representation is limited to the classical and quasi-static 
appr oximat ions . 
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III. FREQUENCY-DEPENDENT RESISTIVITY 

Nevertheless, let us assume that the representation (J2J) can be formally applied also in the 
quantum case. Let us also assume that some prohibition exists for the integrals of the cor- 
relation functions (in the case under consideration the integral of type (jBJ)) to be dependent 
on the " dissipative" parameters (although, as was mentioned above, this requirement is not 
even clearly formulated). Under this assumption it seems that the arguments against the 
CWT in the quantum case in the approximation considered in P] could be correct, because 
the Rq dependence of the integral / is absent in the classical case while it appears in the 
quantum case. 

However, from our point of view, the arguments against the quantum generalization of 
CWT are not correct, because the approximation used in 

[2] is not general. The static 

approximation for resistivity (and, in general, for the impedance) is not correct for high 
frequencies, which are essential in the integral in Eq. (QJ). The dependence R(u) can lead, 
as we will show below, to a dependence of the integral on the effective frequency of collisions 
if L ^ and C ^ oo, even in the classical case. 



A. The classical case for L = 0, C = oo. 

For the special case L = and C = oo we will show below, that: i) the case of R being 
independent of u is meaningless, because in this limit J 2 ~ 1/L — > oo, whereas it must be 
finite, ii) if we use R(u) we arrive at CWT for arbitrary interaction between the particles 
in both classical and quantum cases, iii) there is no fundamental prohibition for correlation 
functions, e.g. the integral J, to be dependent on the " dissipative parameters" (in the sense 
discussed above), although in the classical case this integral in fact is independent of it. 

Now we turn to the main reasons showing that CWT is correct. It was already noted 
that the impedance in quasi-static form, used above, is valid only for very low frequency 
(in the classical state as well as in the quantum case). In particular, for the frequencies 
tkv/T > 1, where the quantum CWT is applicable, the quasi-static form of the impedance 
cannot be used. 

We will now show that even for the classical case the quasi-static approximation for the 
impedance Z(uj), described by Eq. (j2J) cannot be used and, therefore, the Nyquist formula is 
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only an approximation, which is valid only for the quasi-static case of the small frequencies. 
At the same time the CWT is applicable for all frequencies and in this sense is a rigorous 
result for classical and quantum systems. 

For this purpose we consider the specific circuit with L = and C = oo (capacitor with 
the distance between the plates d — 0). Then, in the quasi-static approximation used in Q], 
Z is equal to R in this case and the integral is infinite if R is cj-independent. 

In reality, for this type of circuit (the conductor with L = and C = oo) we have: 

— = ^ (10) 

where cr(u) is the dynamic conductivity, I and S are the length and the area of the conductor 
section, respectively. Eq. (jTUj) is the consequence of the relations j(uS) = <j(u)E(u) and 
J = Sj, e(uS) = lE(uS) for the linear conductor. 

For the simplest case under consideration the approximate expression for the conductivity 
can be written in the Drude-Lorentz form: 

ine 2 

<tdl(w) = —, — — -r, (11) 
m\u) + iv) 

where v is the effective, u- independent, collision frequency of the particles, which transfer 
the current. Substituting Eqs. (|ID jL (jTT )) in Eq. (JJJ) with L = and C = oo we arrive at the 
integral: 

J-oo ^ J-oo 27T 

For the Drude-Lorentz conductivity we find from Eq. (J 12}) 

j2 = 2 DT — Rea DL (u) = / — - T — -. 13 

The result of the integration is, naturally, finite and equals DTne 2 /m = DTQ^/Att, where 
D = S/l. Apparently, this result does not involve the dissipative parameter v. This implies 
that J 2 does not involve the interaction parameter V (we can consider this fact as a specific 
condition for the non-dissipative character of correlation functions, see above). 

It is necessary to stress, that in Eq. (fT2*|l we used the so-called internal conductivity (JUJ), 
which we choose in the Drude-Lorentz approximation. This conductivity is the response 
function for the self-consistent electrical field in Coulomb systems. The general relation 
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between the internal conductivity <r(w) and the external conductivity a ex {u) (the response 
on an external field applied to the system) has the form j^|: 

a{uj) = , . ' , , (14) 

ID 

It is easy to show by straight calculation that the integral ()12|) has the same value for 
and o ex iiS) if we use the Drude-Lorentz approximation. In the general case for the exact 
(non-approximative) expressions for conductivities it is also true Q|, as can be shown by 
using of the Kramers-Kronig relations for a ex [iS) and o[uS). The problem of the validity of 
the Kramers-Kronig relations for conductivities has been investigated by Kirzhnits in [11]. 
In the long- wavelength limit not only a ex , but also the dielectric function e(u) and therefore 
cr(u) satisfy the Kramers-Kronig relations (w.r.t. the external conductivity this statement 
is valid also for finite values of the wave vector k) . The resistivity, present in the relations 
(J2J) , (H21 5 is the internal one, because only the internal conductivity o has a finite limit for 

— ^ and provides a finite value of R'(u —>■ 0) = ReR{uj — > 0). The external coductivity 
in the limit u> — > tends to zero as uj 2 (see, e.g., jl^|). 

By applying the sum rule for the external conductivity of a system of charges with 
arbitrary strong interaction between the particles: 

— Rea ex (u;) = > 15 
-oo 7T ^ m r 

with r the index of the species of particles, we arrive at a general relation similar to Eq. (|12j). 
but without any assumption about the concrete form of the conductivity for the classical 
system: 

J 2 = DT [°° ^Rea ex (uj)=DTY^. (16) 
J-oo 7T *f m r 

The additional assumption that has been adapted to establish Eq. (|15J1 is isotropy of the 
system under consideration. This means, in particular, that the conductivity tensor is diag- 
onal with equal components. According to the statements considered above, Eqs. (fTo^ . (fTp^ 
are also satisfied if instead of a ex we write a. 

As is evident from the general point of view, the current- current correlation function in 
the classical case is independent of the " dissipative" parameters. 



8 



B. The case of a classical system with L / 0, C / oo. 



Let us now consider, as an example of the classical case, the integral © with L ^ and 
C 7^ oo. We assume that the terms related with the inductance and capacitor in Eq. (j2J 
for the impedance are unchanged. For the Drude-Lorentz internal resistivity with constant 
v(u>) = v the real and imaginary parts of the resistivity are equal to ReRi n (u) = R = 
mu/Dne 2 and ImR in {uj) = —muj/Dne 2 , respectively, where ReRi n (uj) is independent of 
oo. For the external resistivity in Drude-Lorentz approximation we obtain ReR(u) = R Q = 
mu/Dne 2 and ImR{u) = -m(uo 2 - Q 2 )/Dume 2 = -Rq(uo 2 - Q 2 )/uuj. 

Then we return to Eq. (JZJ) and to the equality LJ 2 = TL/L with the effective L = 
L — (m/Dne 2 ). It is evident that the condition L > must be fulfilled, which demonstrates 
the limited applicability of the above formulated assumption. However, if we take into 
account that the expression (jl5j) is the result for the thermodynamic limit, whereas the 
Nyquist formula is written for finite samples, we have to consider large values of D and on 
this basis we can suggest L > 0. 

In general, we have to stress the difference between FD relations for finite systems and 
for systems in the thermodynamic limit. Consideration of the case L / and C ^ oo has 
conventional character, because the value L can change sign, although the value J 2 must be 
positive according to definition. This means that representation (JH| has, strictly speaking, 
to be additionally modified for the case L ^ and C ^ oo, when we take into account the 
^-dependence for the resistivity R. 

Nevertheless, it is easy to see, even for the restricted modification of Eq. (|S|) . when we 
use the Drude-Lorentz conductivity with a;-dependent v(uj) the resistivity Rq Rq9{uj) (in 
general, 9{uS) has real and imaginary parts 0'(uj) and 9"{uj)) and the statement about the 
energy LJ 2 being independent of dissipative parameter is not valid. To show this, let us 
consider the conductivity of the Lorentz classical plasma system: 



where fo(p) is the Maxwellian distribution, normalized to the electron density n e , and v(v) = 
AiiZe^neL/rr^v 3 is the velocity-dependent effective collision frequency for charged particles. 
In this approach we use the notation Rq = muo/Dne 2 . The functions 9'(u) and 9"(u) can 
be found by using the relation R = 1 /Da, where a is determined by (|17jl . As follows on 




(17) 



9 



the basis of Eq. ()17j) . the function 9' ^ const and 9" is a nonlinear function of oj/vq. In 
particular, in the high-frequency region (a; > u Q ) the functions 9'{oj) and 9"{u) are equal to: 

32 to 2 
3n uo 2 + v\ 



W = --TT^> (18) 



ou 3 



where v\ = 32z/ /37r. This means that the integral © cannot be represented in a form similar 
to (JTj) with some effective constants L and C and depends on the "dissipative" parameter 

Therefore in the case of a system with L / and C ^ oo the statement about indepen- 
dence of the correlation function of "dissipative parameters" is not true even in the classical 
case, at least for a consideration based on the same dependence Z(u) on L and C as in 
Eq. ©. 

C. The quantum system with L = and C = oo, Callen-Welton relation. 

We now turn to the quantum case with L = and C = oo. The sum rule (II 5 j) is valid 
in the quantum case as well Il3|. At the same time the correlation function J 2 has to be 
written (see, e.g. [131] ) as: 

J 2 = D [ — T w Rea(cj). (20) 

i-oo 

This result can be easily obtained from the general exact CWT for the density current 
correlation function: 

< j,(r, t)j k (r', t') >= ^ coth f^j [a ik (uj) + a* ki (u;)) S(r - r') (21) 

where o~ik is a component of the tensorial conductivity and the asterisk denotes the complex 
conjugate by integration over r and r' in the volume of a linear conductor V — IS and the 
suggestion that in this volume the conductivity is the same as in the thermodynamic limit. 
Under this suggestion the result coincides with the Nyquist relation in the form (|2(Jj) . where 
conductivity is an exact function of the frequency to. 

There exist no arguments for independence of this integral of the parameter V. Therefore 
this integral depends on "dissipative parameters". Let us demonstrate this fact in the 
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example of the expansion of J 2 in the parameter Tiuo/T. In this case the Eq. ()20|) can be 
rewritten approximately in terms of the zeroth and second momenta of the conductivity: 



oc 



J^Dj ^T.ReaM^^o + ^), (22) 

7T 7T 111 A 



where 



/oo 
duju 2n Rea(uj). (23) 
-oo 

The term with fi = leads to the classical value of J 2 (see Eq. (J20|) ). The second 

moment of the external conductivity has been calculated in 10| and can be expressed in 
terms of the value of the electron- ion equilibrium correlation function g e i(r) at r = 0. It can 
be connected with the second moment of the internal conductivity by use of the asymptotic 
expansions for the Kramers- Kronig formulas. For the case Z = 1 the second moment of the 
internal conductivity equals: 

The final expression for J 2 following from the relations mentioned above is: 

Naturally the function g e i{r), as well as it value in the point r = 0, is a function of the 
parameter T. To calculate g e i(0) for a purely Coulomb system the interaction and quantum 
effects in g ei also have to be taken into account, in particular to avoid the divergence g e i(0) 



for small distances. For example, 
function can be written (Z=l) as 



or a weak plasma interaction T 1 the pair correlation 



3 



9ei(r) = S£ 2) (r)exp 



,2 



(26) 



AirTr 

where k^ 1 = ^JT /Aiie 2 n is the Debye radius (with the density n = n e + Hi) and S^\r) is 
the pair electron-ion Slater sum, which takes into account the quantum effects. The value 
g e i(0) for the approximation (J20"|) )equals: 

/ F3/2 \ 

^(0) = 5i 2) (0)exp(^- 7 =J, (27) 
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The Slater sum Sj? (0) is a function of the parameter ( = e 2 / 4irT \ ei , where A e j = h/ \JlmT 
and m are De Broglie wavelength and the reduced electro-ion mass respectively. This func- 
tion has a finite value for all values of the parameter X ei , except A e j = 0. 

Therefore J 2 , even in the simple approximations used above, depends on the " dissipative" 
parameter T. The dependence of the current- current equilibrium correlation functions on the 
"dissipative" parameters is not in contradiction with any fundamental principles of quantum 
statistical theory. 



IV. CONCLUSIONS. 



The main conclusion of the present paper is the statement that the quasi-static approx- 
imation for the impedance in the form used in ^, |^ is not allowed when considering the 
opportunity for violation of the general quantum CWT. Therefore, presently there are no 
arguments to discuss violation of the CWT in the quantum case on this basis. As it is, 
strictly speaking, the CWT is proved only for the case of a conductor with L = and 
C = oo and for an infinite medium, when the Hamiltonian description of the system is well- 
determined. For this case we have shown the way to connect the macroscopical model of 
the Nyquist fluctuation-dissipation relation with the microscopic derivation of the CWT by 
introduction of the frequency dependent resistivity in the Nyquist model. We also demon- 
strated the absence at this time of a clear way for generalization of this relation for the 
case L ^ and C ^ oo, in spite of some attempts to do that ^j]. The important question 
about determination of the "dissipative" parameters has been discussed above for plasma- 
like systems and connected, in particular cases considered above, with the T-dependence of 
the current- current correlation function. In general, a more fundamental determination of 
the dissipative parameters has to be formulated. 

In the framework of the kinetic equations, the dissipative parameters are connected with 
the respective effective collision frequencies. For the problem of FD relations the approach 
of kinetic equations is not enough, because the results have to be found for arbitrarily strong 
interaction between the particles. In this case the determination of the dissipative parameters 
is based, in general, on the imaginary parts of the poles of the analytical continuation of the 
temperature Green functions In the particular cases of weak or short-range interactions 
between the particles these poles lead to the same results as the kinetic equations approach. 
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The problem of the validity of summation of the rows of the perturbation theory for the 
Green functions, as well as the problem of irreversible behavior of physical systems due 
to damping, connected with these poles, have a fundamental character and still have not 
been completely clarified. Nevertheless, we trust that the general results, connected with 
dissipation based on the Green function approach are applicable to obtain the fluctuation- 
dissipation relations for systems with well-defined Hamiltonian. Of cause the T- dependence 
of the current-current correlation function, which we consider as a manifestation of the 
dependence on "dissipative parameters", can be related in fact not only with a shifting of 
the imaginary part of the poles of the Green functions, but also with a shifting of the real 
parts of these poles. Nevertheless we do not see the opportunity to distinguish in the exact 
final result for CWT the "dissipative" and "non-dissipative" dependence on the parameter 

r. 

The sense of the consideration, suggested by Yu.L. Klimontovich \2\ is, from our point of 
view, not about incorrectness of the CWT, but about limited applicability of the Nyquist 
formula in the approximation of constant R. This limitation follows, at first, as was shown 
above, from the quasi-static character of the impedance in the form (j2J). It is evident not 
only in the quantum, but also in the classical case. At the moment there is no clear basis for 
a general and rigorous consideration of the electrical circuit with finite L ^ and C 7^ 00 
in the quantum and even in the classical case, except the model of quasi-static approxi- 
mation for the impedance, which is applicable nevertheless to many experimental situations. 



can be calculated analytically. According to the Maple calculation program this integral 
equals one. We can rewrite the integral (|28jl in the form 



V. APPENDIX. 



The integral I in Eq. (JZj) 




(28) 




(29) 
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with A = L/CR^ > 0. The denominator is equal to (£ 2 — £ 2 )(£ 2 — £f), where the roots are: 



£ 2 = - [-(1 - 2A) + ^/il^AX) 



5 S2 



-(1-2A)- V(1-4A) 



The integral ((Zj) can be split into two integrals I = 1\ + 1%, with 



h 



di- 



(30) 



31) 



TTv^r^Ayy-oo "(e-e?)' tt vo^^a) ./-oo ^e-e 2 2 ) 

Let us consider the particular case 1 — 4A > 0. It is easy to see that the roots £ 2 ,£f in this 
case are real and negative. Therefore the integrals are equal to: 



h 



(32) 



V^4A)' y^T^X)' 
Because | ^| |>| ^ 2 | the sum I = Ii + I 2 > and equals unity. To show this it suffices to 
take I 2 and to use Eq. (132(1 . The same is true for an arbitrary value of A. 



Acknowledgment 



Authors are gateful to W. Ebeling, A.M. Ignatov, N.I. Klyuchnikov, and I.M. Tkachenko 
for valuable discussions of the various problems, reflected in this work, which has been 
submitted in Physica A (8 March 2006). 



[1] L.D. Landau, E.M.Lifshitz, Statistical Physics, "Science" Publishing, Moscow, 1976. 
[2] Yu.L. Klimontovich, Uspekhi Fiz. Nauk, 151, 2, 309 (1987). Private communications (1999- 
2001). 

[3] B.A. Veklenko, Yu.B. Scherkunov, Applied Physics (in Russian), 4, 5 (2004). 

[4] H. Nyquist, Phys. Rev., 32, 110 (1928). 

[5] G.S. Gorelik, Uspekhi Fiz. Nauk, XLIV, 1, 33 (1951). 

[6] H.B. Callen, Th.A. Welton, Phys. Rev., 83, 34 (1951). 

[7] L.D. Landau, E.M. Lifshitz, Electrodynamics of the Continuous Media, Fizmatgiz, Moscow, 
1959. 

[8] S.A. Trigger, Meetings with Yu.L. Klimontovich, Condensed Matter Physics, 7, N 3 (39), 443 
(2004); DAMU-Hefte, 1 (2005) 56 (Berlin). 



14 



[9] A.I.Akhiezer, S.V. Peletminski, Methods of Statistical Physics, "Science" Publishing, Moscow, 
1977. 

[10] V.M. Adamyan, I.M. Tkachenko, Contr. PI. Phys. 43, 252 (2003). 
[11] D.A. Kirzhnits, Uspekhi Fiz. Nauk, 119, 2, 357 (1976). 

[12] D. Pines, P. Nozieres, The Theory of Quantum Liquids, v.I, W.A. Benjamin, INC. New York 

- Amsterdam, 1966. 
[13] R. Kubo, J. Phys. Soc. Japan, 12, 570 (1957). 

[14] W.D. Kraeft, D. Kremp, W.Ebeling, G. Roepke, Quantum Statistics of Charged Particle 

Systems, Akademie Verlag, Berlin, 1986. 
[15] A. A. Abrikosov, L.P. Gor'kov, I.E. Dzyaloshinskii, Methods of Quantum Field Theory in 

Statistical Physics, Fiz-Mat. Publishing, Moscow, 1962. 



